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Energy-dependent orbital modulation of X-rays and constraints on 
emission of the jet in Cyg X-3 
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ABSTRACT 

We study orbital modulation of X-rays from Cyg X-3, using data from Swift, INTEGRAL and 
RXTE. Using the wealth of the presently available data and an improved averaging method, 
we obtain energy-dependent folded and averaged light curves with unprecedented accuracy. 
We find that above ~5 keV, the modulation depth decreases with the increasing energy, which 
is consistent with the modulation being caused by both bound-free absorption and Compton 
scattering in the stellar wind of the donor, with minima corresponding to the highest optical 
depth, which occurs around the superior conjunction. We find a decrease of the depth below 
~3 ke V, which appears to be due to re-emission of the absorbed continuum by the wind in soft 
X-ray lines. Based on the shape of the folded light curves, any X-ray contribution from the jet 
in Cyg X-3, which emits y-rays detected at energies > 0.1 GeV in soft spectral states, is found 
to be minor up to ~100 keV. This implies the presence of a rather sharp low-energy break in 
the jet Me V-range spectrum. We also calculate phase-resolved RXTE X-ray spectra, and show 
the difference between the spectra corresponding to phases around the superior and inferior 
conjunctions can indeed be accounted for by a combined effect of bound-free absorption in 
an ionized medium and Compton scattering. 
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1 INTRODUCTION 

Cyg X -3 is a high-mass X-ray binary with a Wolf-Rayet (WR) 
donor ( Ivan Kerkwiik et alJI 19921, 1 19961 : Ivan Kerkwiildl 19931) and a 
very short orbital period of P ^ 4.8 h. I t s distance is ^ 7-9 kpc 
dLing. Zhang & Tandl2O09l ; iDickeyll 19831: Predehl et ai]|200d) . In 



spite of its discovery already in 1966 dGiacconi et al 119671) . CveX- 
3 remains poorly understood. In particular, the nature of its compact 
object remains uncertain, due to the lack of relia ble determination 
of the mass functions and inclination (see, e.g.. IVilhu et al.ll2009l 
for a discussion). However, the presence of a black hole is favoured 
by consider ing the X-ray and radio emission and the bolometric 
luminosity (Hjalmarsdotter et al. 2008, 2009; Szostek & Zdziarski 
l2008L hereafter SZ08, ISzostek. Zdziarski & McColloughl I2008L 
hereafter SZM08). Also. IZdziarski. Misra & GierlinskT d2010l) 



have shown that the differences in the form of the X-ray spectra 
of Cyg X-3 from those of confirmed black-hole binaries can be ac- 
counted for by Compton scattering in a cloud formed by the stellar 
wind from the companion. 

Cyg X-3 is a persistent X-ray source with a typical X-ray lu- 
minosity of L x ~ 10 38 erg s _I . Its X-ray spectra have been classified 
into five states by SZM08, who have also quantified their correla- 
tions with the radio emission. Its high-energy y-ray emission has 
been discovered by the Fermi Lar ge Area Telescope (LAT) and by 
AGILE in the sof t spectral states ( Fermi LAT Co llaborationl l2009L 



hereafter FLC09; iTavani et~aT 



and A GIL E are presented b y 



20091). Later d e tectio n s by the LAT 



Willia ms et ail (1201 ll) . ICorbel et al.l 



d2012l) and lBulgarelli et al.ld2012l) . hereafter B 12. 

Cyg X-3 shows pronounced flu x modulation o n 

4.8 h period, discovered in X-rays dParsignault et all 1 19721; 



the 



ISanford & Hawkins! 1 19721; ICanizares et al.l 1 1973b and 
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frared dBecklin et al.l 19731) . Both periodicities are most likely 
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related to the strong stellar wind from the WR star (e.g., 
IWillingale. King & Poundsll 19851 ; Ivan Kerkwiidl 19931) . In particu- 
lar, the X-ray flux minima correspond to the maximum absorption 
in the line of sight around the superior conjunction (the compact 
object be hind the WR star). The X-ray period has been found to 
increa s e dManzo. Molteni & Robballl978l; lLamb. Dower & Fickle! 



1979; Eisner et al 



19801; Ivan der Kris & Bonnet-Bidaud Il98ll . 



19891 ; IKitamoto et alj|l98lll995l ; Isingh et alj|2002l) . The period 
increase of PjP m 1 x 10~ 6 yr~' is likely to be due to the loss of an- 
gular momentum through the wind mass loss jPavidsen & Ostrikej 
fl974h . 

In this work, we present a detailed study of the dependence 
of the form of the X-ray orbital modulation on energy, using data 
from Swift, INTEGRAL and Rossi X-ray Timing Explorer (RXTE). 
We also study the dependence of the modulation on the X-ray spec- 
tral state, as well as obtaining X-ray modulation profiles for inter- 
vals simultaneous with the detections of high-energy y-ray s. The 
y-ray emission was found to be also strongly modulated on the or- 
bital period, but with the maxima and minima approximately corre- 
sponding to the X-ray minima and maxima, respectively (FLC09). 
This phase reversal has been interpreted as due to anisotropy of 
Compton scattering of the stellar emiss ion by relativistic electrons 
in the jet jDubus. Cerutti & H enri 2010). The jet emiss ion certainly 
extend s to lower energies, and, using the model of Dub us et al.l 
(2010), the hard X-rays from the jet are also predicted to have the 
maximum around the superior conjunction dZdziarski et al.ll2~012l 
hereafter Z 12), i.e., shifted with respect to the observed X-ray mod- 
ulation by a half of the period. A measurement of the energy- 
dependent orbital modulation up to hard X-rays simultaneous with 
the y-ray emission can then constrain the X-ray contribution of the 
jet. 



2 FOLDED LIGHT CURVES FROM X-RAY 
MONITORING 

2.1 The X-ray monitoring data 

We use X-ra y monitoring data from the Swift Burst Alert T ele- 
scope (BAT; iBarthelmv elat] 120051 ; iMarkwardt et al.l l2005h in 
the form of a 14-195 keV 8-channel light curve created for 
this study. The typical exposure for a single data point is ~ 
10 3 s. The data contain 23783 points within 1807 days span- 
ning MJD 53355-55469. The channels are between energies of 
14, 20, 24, 35, 50, 75, 100, 150 and 195 keV. These data 
have been obtained using an analysis different from that of 
http://swift.gsfc.nasa.gov/docs/swift/results/transients , which gives 
BAT count rates in the single 15-50 keV channel. Fig. Q] shows 
a comparison of the daily-averaged 15-50 keV count rates with 
the corresponding 14-50 keV energy fluxes from our analy- 
sis. The energy fluxes have been obtained by converting the 
count rates of our data using sca ling to the Crab spectrum (see 
Izdziarski. Poolev~& Skinner 2011). 

We also use both dwe ll and daily-averaged data from 
the All-Sky Mon itor (ASM; iBradt, Rothschild & Swankl 1 19931 : 
iLevine et alll 19961) on board RXTE. The dwell data contain 97996 
measurements within 5267 days spanning MJD 50087-55915, and 
the exposure of a single observation is < 90 s. The ASM has three 
channels at energies of 1.5-3 keV, 3-5 keV and 5-12 keV, for which 
the corresponding energy fluxes are obtained by scaling to the Crab 
spectrum, as above. 

Fig. [2] shows the relationship between the energy fluxes in 
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Figure 1. The correlation between the daily averages of the BAT flux, F(14- 
50 keV) from our analysis vs. those of the 15-50 keV rate in the public BAT 
data, see Section |2~T| Only points with statistical significance > 3o" (for each 
of the quantities) are shown. The black cross shows the flux/rate of F(14-50 
keV) = 1.2 keV cm" 2 s" 1 , corresponding to R(15-50 keV) = 0.028 cm" 2 
s , below which we define our soft spectral state. 
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Figure 2. The correlation between the daily-averaged energy fluxes in the 
3-5 keV and 14-50 keV ranges. The hard and soft states are defined here by 
F(3-5 keV) < 0.58 keV cm" 2 s" 1 and F(14-50 keV) < 1.2 keV cm" 2 s _I , 
respectively. The boundaries are marked by the dotted and dashed lines. An 
intermediate state corresponds to both the 3-5 keV and 14-50 keV fluxes 
above the respective boundaries. Only points with statistical significance 
> 3cr (for each of the fluxes) are shown. 
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the 3-5 keV and 14-50 keV photon energy ranges. We see that 
at low soft X-ray fluxes, about F(3-5 keV) < 0.58 keV cirT 2 s _1 
(marked by the dotted line), corresponding to the ASM 3-5 keV 
count rate < 3 s~', there is a clear soft/hard X-ray flux anticorre- 
lation. Note that a positive correlation between the 3-5 keV count 
rate and radio flux is seen below the same rate (SZM08). At higher 
3-5 keV fluxes but at high 14-50 keV fluxes (F( 14-50 keV) > 1.2 
keV cirT 2 s" 1 , marked by the dashed line), the hard X-ray flux 
changes within a narrow range independent of the soft X-ray flux. 
At lower 14-50 keV fluxes (below the dashed line), there is an ap- 
parent anticorrelation down to F(14-50 keV) ~ 0.5 keV cirT 2 s~'. 
Below it, the soft X-ray flux changes within a narrow range approx- 
imately independent of the hard X-ray flux. Such anticorrelations, 
expressing a spectral pivoting around ~ 10-15 keV, can also be in- 
ferred from figs. 1-2 of SZ M08. Fig. 7b of SZM08 and fig. 3 of 
Hialmarsd otter et all J2008) show similar correlations with the 20- 
100 keV flux measured (though with a lower sensitivity than that of 
the BAT) by the Burst and Transient Source Experiment (BATSE) 
on board Compton Gamma Ray Obse rvatory (CGRO) . 

We have also studied MAXI dMatsuoka et alj 2009) data. 
However, the exposure of a single data point for those data is 1.5 h, 
which is ~ P/3 and much longer than any orbital bin we use. Thus, 
we do not use those data for studying orbital modulation in Cyg X- 
3. We note those data also provide us with a hard-state condition, of 
the 2^1 keV photon flux < 0.2 cirT 2 s~', approximately equivalent 
to the ASM 3-5 keV count rate < 3 s -1 . 

We performed analyses including all of the ASM data and sim- 
ilarly with all of the BAT data, as well as some subsets of each. 
First, we use a hard-state criterion of the daily-averaged 3-5 keV 
ASM count rate < 3.0 s~', corresponding to F(3-5 keV) < 0.58 
keV cm" 2 s _1 , to the left of the dotted line in Fig. [2] As the soft- 
state criterion, we use the daily-averaged 14-50 keV BAT flux of 
< 1.2 keV cirT 2 s -1 , i.e., below the dashed line in Fig. [2] Our inter- 
mediate state corresponds to both fluxes being above these bound- 
aries, i.e., in the upper right corner of Fig. |2] The soft-state cri- 
terion corresponds to the intervals of MJD 53746-53896, 53913- 
54138, 54204-54252, 54445-54470, 54548-54642, 54748-54818, 
54980-55044, 55080-55103, 55322-55349 and 55584-55645 (de- 
termined using a 5-d running average in order to reduce the effect 
of fluctuations). For epochs after the end day of our BAT data, we 
have applied the corresponding criterion for the public 15-50 keV 
BAT data. 

We also use data from intervals corresponding to the detected 
high-energy y-ray emission. Currently, they are MJD 54566-54647 
(B12), 54750-54820 (FLC09), 54821-54850 (B 12) 54990-55045 
(FLC 09), 55324-55326 (B12), 55343-553 45 ^Williams et al.1 
l201ll ; B12) and 55586-55610, 55642-55644 dCorbel et alj|2012l) . 
All those periods correspond to the soft state as defined above, i.e., 
below the dashed line in Fig. [2] see FLC09 and B12. Thus, the ap- 
pearance of y-ray emission corresponds to a marked change of the 
form of the soft/hard X-ray correlation. Combining the results from 
Fermi LAT and AGILE, it appears likely that all low hard-flux in- 
tervals are also associated with some high-energy y-ray emission. 



2.2 Treatment of the light curves 

Cyg X-3 has an orbital period that is increasing with time. Thus, 
we use a quadratic ephemeris, 



time, To, Po is the period at To, P is the peri od derivative, and P is 
the period at T,„. The ephemeris (in UTC) of Si ngh et all d2002l) is 

T = 40949.392 [MJD], P = 0.19968443 d, c = 5.75 x 10 -11 d.(2) 

An updated ephemeris taking into account Suzaku observations (S. 
Kitamoto, private communication) is 

T = 40949.3913 [MJD], P = 0.19968451 d, c = 5.62x10"" d,(3) 

in which case P/P =s (1.03 ±0.02) x 10" 6 yr 1 . 

We note that these ephemerides of Cyg X-3 use the template 
of lvan der Klis &B onnet-Bidaud fl989» . which is defined numeri- 
cally by their table 2. That template has the minimum at the phase 
- 0.96-0.97 instead of 1.0 (or, equivalently, phase 0.0). Thus, min- 
ima at phases < 1 obtained by us do not indicate a discrepancy 
with respect to previous results. We also correct the light curves for 
barycentric delays. 

We need to properly average the folded light curves, [F,(£;), 
o~i(ti)], where F, is the average count rate of the observation with 
the mid-time of r,, and cr, is its measurement error. One issue to 
consider is the large dynamic range of the flux variability of Cyg 
X-3, which is almost two orders of magnitude (see Fig. |2}- Thus, 
linearly-averaged modulation profiles are strongly biased by those 
at highest flux states. Using flux logarithms can, in principle, some- 
what alleviate this problem, as it strongly reduces the dynamic 
range used in summation. However, this requires removal of all 
negative flux measurements. For a weak signal, this introduces a 
strong bias against the minima of the modulation. We found that 
this effect strongly suppresses the depth of the modulation at ener- 
gies > 30 keV. 

On the other hand, orbital modulation acts on a signal with 
a local flux level, the level being governed by the source aperiodic 
variability. The periodic and aperiodic variabilities are largely inde- 
pendent of each other. Thus, in order to effectively reduce the large 
amplitude of the underlying aperiodic flux variability, we calculate 
a running average, averaging all observations within some ±At of a 
given observation time, 



< Af, 



(4) 



To + Pom + cm 1 , c = P P/2, P = P + 2cm, 



(1) 



where T m is the time of a ra-th occurrence of a zero orbital phase 
(presumably related to the superior conjunction) since the reference 



where J is the number of data points satisfying the above condi- 
tion on j. Then, we renormalize both F, and cr, by dividing each by 
F (so CTj/F, remain unchanged). In this way, effects of the aperi- 
odic long-term variability on time scales > 2At are removed. In the 
case of measurements equally spaced at a time interval of <k Pj K, 
where K is the number of bins per period at which the folded light 
curve is calculated, using At = P/2 would optimally remove the ef- 
fect of long-term aperiodic variability. However, our measurements 
do not satisfy this criterion, and are occasionally sparsely and un- 
evenly spaced. Then, the shorter the At, the worse the statistics on 
which the running average is based. In the extreme case of no other 
measurements within ±At from a given observation time, F, would 
be reset to unity, and such a point would falsely reduce the actual 
orbital modulation. Thus, in our method, we reject points for which 
the running average is based on J < 10 points. We have tested this 
method for our data, and have found that At = 1 d (^ 5F) is ap- 
proximately the shortest interval for which the running average is 
estimated with good statistics and only a small fraction of points 
have J < 10. Namely, (J) ± rms, and the fraction of the rejected 
points are (51 ± 23, 0.02), (34 ± 16, 0.04), for the ASM and BAT 
data, respectively. We thus use the above value of At hereafter. 

Another issue here is related to aperiodic short time-scale vari- 
ability on time scales shorter than the length of an orbital bin, which 
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effect we also would like to remove. Therefore, we pre-average the 
light curves on real-time intervals with the length equal to the bin 
size, P/K, where K is the number of phase bins per orbit. We obtain 
local averages of the renormalized (see above) flux and its error, 
which we denote as F^, a>;, respectively, where k is the number 
of the phase bin and / is the number of the time bin (counted from 
the start time of a given data set) contributing to the k-th phase bin. 
This m ethod was applied in studying pe r iodic variability of Cyg 
X-l by Ibragimov, Zdziarski & Poutanen (2007). It also partly re- 
moves a bias on folded averages due to a non-uniform coverage, 
e.g., due to a large number of points during one time interval and a 
low number of points during another, with both contributing to the 
same phase bin, k. Without pre-averaging, the former would have 
much higher weight than the latter, even if all measurement were 
relatively accurate. The pre-averaging does not have a major effect 
in the case of Cyg X-3 due to its very short orbital period (result- 
ing in rather few points per a P/K time interval), as well as for 
ASM and BAT data, which have relatively uniform coverage. Still, 
this is the statistically correct procedure, resulting, in particular, in 
correctly estimated uncertainties of the folded and averaged light 
curves. 

In final averaging, we use weights given by inverse squares of 
the uncertainties, 



F k ,lai, 



T Nl l/o- 2 



(5) 



where Nu is the number of time bins contributing to the fc-th phase 
bin. The error of the weighted average is estimated from the vari- 
ance, {{F 2 k ) - {F k ) 2 )/(N k - 1), 



y Nk II a- 2 



' (F k Y 



(6) 



We fit the obtained folded/averaged light curves by a sim- 
ple model of absorption/scattering in a spherically symmetric stel- 
lar wind. A commonly-used wind velocity profil e from a massive 
star (e.g. jLamers. Cerruti-Sola & Perinottoll 987) and the resulting 
electron density are, respectively, 



-M 



(7) 



4nm v fx e r 2 v(r) ' 

where r is the distance from the centre of the donor of the radius R„ 
Voo is the terminal velocity, f} parametrizes the wind acceleration, M 
is the total mass-loss rate, m p is the proton mass, yu e 2/(1 + X) 
is the mean electron molecular weight, and X is the H fraction. A 
small correction due to velocity reaching the sound speed rather 
than being null at the stellar surface has been neglected. 

Assuming the opacity for a given photon energy is a spatially 
independent constant times the density (as is the case, e.g., for 
Compton scattering or for bound-free absorption with a constant 
ionization c oefficient), th e optical depth is an integral over the pho- 
ton path, / dPringldll974l) . 

r(0) = r o f (_) 2 fi-^) da/a), (8) 



where tq = (a"c + cr mn )noa, a~c is the Compton cross section (in- 
cluding its Klein-Nishina decline), cr ion is a photoionization cross 
section per electron, n is a fiducial electron density at r = a un- 
der the assumption of v = Vo and a is the orbital separation. For a 
circular orbit and assuming the radiation originating at the compact 
object location, r is related to I by 



where i is the binary inclination and <f> is the orbital phase (defined 
as (p = 0, n at the superior a nd inferior conjunction, respectively), 
see fig. 1 in Zdziarski fe012h . For r ^> 7?*, i.e., v — v M , 



r(0) - t 



tt/2 + arcsin s 



(10) 



which is is equivalent to equation (8) o f |Pringldll974T) . 

In Compton scattering, a photon removed from one line of 
sight appears in a different one. Photons scattered away from di- 
rections with high r preferentially leave the system close to the 
directions with low t. Also, the photon distribution will be a func- 
tion of the angle with respect to the normal to the binary plane. It 
is possible to model this effect quantitatively in a more detailed 
treatme nt, e.g., using a Monte Carlo m ethod. Such method was 
used bv lHertz. Joss & Rappaporq dl978h . who, however, presented 
their results for i = 90° only (and for isotropic wind). However, 
the wind density distribution in Cyg X-3 is certainly not isotropic, 
but instead with a str ong focusing towards the compact object 
dFriend&Castonll982l) . and the wind structure is affected by ir- 
radiation from the X-ray source. Taking into account these compli- 
cations is beyond the scop e of this work . We therefore make a sim- 
plifying assumption (as infPringle 1974]) that a scattering removes 
the photon from the line of sight and neglect photons scattered into 
it. This allows us to fit the observed phase-dependent fluxes by 



F((P)/F = exp[-T(0-0 o )], 



(11) 



with To, i, j8, the offset phase, 0o, and Fq being the model free pa- 
rameters. We note that F = F max exp[r(7r)]. This requires the scat- 
tering optical depth to be < 1 (though the absorption optical depth 
can be any). We define the fractional modulation depth as 



D = 



1 - exp(-Ar), At = t(0) - t(7t), 



(12) 



(r/a) = 1 + (I I a) - 2(//a)sin;cos0, 



(9) 



where F max and F m ; n are the maximum and minimum modula- 
tion fluxes, respectively. For the model of equation dlOt . At ^ 

2Toi/ cos i. 

We then need to specify the parameters of Cyg X-3 appear- 
ing in this model. Since Cyg X-3 is a He star, X ^ and yu e =^ 2. 
T he stellar radius, /?», i s clos ely related to the WR mass, see fig. 2 
of ISchaerer & Maeded dl992l) . The separation, a follows from the 
total mass of the system. This, however, remains rather uncertain. 
Here, we choose, so mewhat arbitrarily, a binary solution satisfying 
all the constraints of IVilhu et al] (2009), namely M, = 20Mq and 
the compact object mass of Mc = 6Mq, for which a ^ 3.0 x 10" 
cm. T his M, corresponds to R, ^ 1.0x10" cm dSchaerer & Maederl 
1992). This solution corresponds to an allow ed range of i nclina - 
tions of about i =s (33°-63°), see section 6 of IVilhu et al] d2009h - 
A typical value of the terminal wind velocity estimated for Cyg 
X-3 is v„ - 1.7 x 10 s cm s~' (see a discussion and references in 
SZ08). We assume fi = 2, which was fou nd to fit well the re sults of 
wind particle simulati ons for Cyg X- 3 of IVilhu et al] d2009l) . Also, 
this value was used bv lLanger! dl989f) in his modelling of mass loss 
from WR stars. 

As an alternative, and purely phenomenological, model for the 
folded light curves, we use a sum of three Fourier harmonics (in the 
logarithmic spac e). We follow here exa ctly the treatment described 
in section 3.2 of lLachowicz et all d2006h . The modulation depth in 
this model is determined numerically. 
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Figure 3. (a) Folded/averaged light curves for the count rates for the energy 
range of 1.5-12 keV (blue, ASM), 14-50 keV (red, BAT) and 20^0 keV 
(green, INTEGRAL/1SGRI), using K = 40, 20 and 12 bins, respectively. For 
the BAT data, only observations with the exposure < 0.5P/K have been se- 
lected. Here and elsewhere, two orbital ranges are shown for clarity of the 
display, (b) Comparison of the BAT folded light curve with the above expo- 
sure criterion and K = 20 with that without it (for K = 40). (c) Comparison 
of the ASM folded light curve calculated with our treatment described in 
Section [2~2l (blue triangles), i.e., using the barycentric-corrected data, the 
quadratic ephemeris of equation (|3j, the renormalization of the light curve 
to the running average, and the pre-averaging, with that obtained using a 
constant orbital period, unweighted averaging and no other corrections (red 
squares). 



2.3 Energy-dependent folded light curves 

We first consider folded/averaged light curves averaged over all the 
available ASM and BAT data. We have compared the light curves 
folded with the ephemerides of equations (0 and ((3}- We have 
found the differences are tiny, and we use the updated ephemeris 
of equation l[3} hereafter. 

Fig. [Ha) shows the folded light curves for the total count rates 
of the ASM and BAT instruments (and also shows the 20-40 keV 
profile from INTEGRAL, see Section [3] below). We see the frac- 



-i — | — i — i — i — i — | — i — i — i — i — | — i— 
ASM all 1.5-3 keV 




0/2n 

-i — | — i — i — i — i — | — i — i — i — i — r 
BAT all 14-20 keV 




0/27T 

Figure 4. Folded light curves using all of the data for (a) 1.5-3 keV ASM, 
and (b) 14-20 keV BAT data, using K = 40. The red solid curves show fits 
to the data with the model of equation {8), (TTJ. The green dashed curves 
show fits by the 3-harmonic model. 



1.4 
1.2 

A 

1 

V 1 

\ 

0.8 
0.6 
0.4 



1.4 

1.2 

A 

V 1 

\ 
&. 

0.8 
0.6 
0.4 



i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 

ASM hard/intermediate/soft 1.5-12 keV 

.t !! ». 

i i 

; i 



0.5 



1.5 



0/2n 



— , , | , , , , | , , , , | , , 
BAT hard/intermediate/soft 14-50 keV 



8 • 



0.5 



1.5 



0/27T 



Figure 5. Comparison of the folded light curves in the hard (blue triangles), 
intermediate (green squares) and soft (red circles) spectral states for (a) the 
ASM 1.5-12 keV range, and (b) the BAT 14-50 keV range. 
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Figure 6. Comparison of the folded light curves in different energy channels 
for all the available ASM and BAT data. The blue, red and green light curves 
are for channels (a) 1 .5-3 keV, 3-5 keV, 5-12 keV, (b) 5-12 keV, 14-20 keV, 
50-75 keV, and (c) 24-35 keV, 75-100 keV, 100-150 keV, respectively. 



tional modulation depth, equation ( 112b . is significantly larger in the 
1.5-12 keV range, D =* 0.61, than the 14-50 keV one, D 0.50. 
This is likely due to the bound-free absorption by the wind being 
stronger in the ASM band than in the BAT one. The BAT minimum 
is at (f>/2n m 0.95, somewhat less than the ASM one at ^ 0.97. 
Both profiles have relatively complex and asymmetric shapes, with 
a flattening around the phase of 0.4 followed by a peak around 0.55- 
0.60. The ASM profile has two statistically significant dips, at the 
phases 0.40 and 0.575. For either of the dips, the difference with 
respect to the rate averaged over the two neighbouring bins is > 5 
standard deviations. Their presence indicates some complexity of 
the wi nd structure, pro bably due to focusing towards the compact 
object jFriend & CastoJ 19821) . 

An issue for the BAT data concerns their relatively long typi- 
cal exposure of a single data point of ~ 10 3 s. This is comparable 
to the length of one orbital bin in the case of K = 20, which may 
cause a smoothing of sharp features in the folded light curves. In 



Fig. [^a), we have thus used only measurements with the exposure 
< 0.5P/K for the BAT profile. However, using only BAT data with 
short exposure times strongly reduces the number of used measure- 
ments. In the present case, accepting only exposures < 0.5P/K for 
K = 20 reduces the number of available measurements by a factor 
of ~ 10 with respect to using all the available data. In Fig.[3{b), we 
compare the profiles with and without imposing this condition. We 
see that using the data without this selection causes certain smooth- 
ing of the profile, but the effect is rather minor. Thus, we do not use 
this selection for the BAT data hereafter. 

Fig.[3{c) illustrates how taking into account various effects can 
affect folded/averaged light curves, for the example of the ASM 
light curve. The red squares show the profile obt ained with a con- 
stant P computed for the ASM data up to 2004 bv lWen et ail ( 2006) 
(i.e., neglecting P), neglecting the barycentric correction, without 
prior renormalizing with respect to the running average and with- 
out pre-averaging within time bins. We see that that curve has the 
minimum at <pl2n ^ 0.9 instead of 0.95 found for the treatment 
of Section I2T2I As we have checked, this shift is almost entirely 
due to neglecting the P of Cyg X-3. Taking into account barycen- 
tric time delays introduce only a very small effect. Furthermore, the 
red points have a muc h flatter and more scattered top part (similar 
to the profile shown in lWen et al.ll2006ri . This difference is mostly 
due to renormalizing the light curve w/r to the running average. The 
profile without it is strongly dominated by the brightest states (as 
pointed out in Section [2~2t . Then, fluctuations of those dominant 
states have a strong effect. On the other hand, our profile, shown by 
the blue triangles, is evenly averaged over all flux states. Finally, 
the pre-averaging introduces only a small effect for the ASM data, 
though it does noticeably increase the error bars. 

The 14-50 keV BAT profile shown in Fig. [3] can be com- 
pared with the folded light curve in the 12-50 keV band from the 
Gam ma-ray Burst Monito r (GBM ) on board Fermi, shown in fig. 
11 of IWilson-Hodge et alj ( 120121) . Unlike our, rather smooth, BAT 
profile, the GBM curve has multiple s harp peaks. We note that the 
orbital profile for Cyg X- 1 shown in Wil son-Hodge et all d2012l) 
also has much stronger sharp fea tures than the corresp onding 15- 
50 keV BAT profile presented by IZdziarski et all d201 ll) . Possibly, 
the sharp GBM features are due to statistical fluctuations at some- 
what underestimated measurement uncertainties. 

In Fig.[4{a-b), we fit the lowest-energy channels of the ASM 
and BAT with the model of equations {8]l, i ll 1 b . We see that neither 
profile is well fitted by this model. The profiles are significantly 
asymmetric, with a slower rise and a faster decline. (The scatter 
around the top of the 1.5-3 keV profile appears statistical.) Thus, 
the model does not provide statistically good fits, but we give the fit 
parameters for completeness: i = 33°±1°, 39°±1°,t = 0.27±0.01, 
0.15 ± 0.01, 4>ol2n = 0.02 ± 0.01, -0.03 ± 0.01, for the ASM and 
BAT, respectively (for the fixed /? = 2, see Section 12.21 ). Given 
the statistically unsatisfactory fits, these parameters may not cor- 
respond to the actual wind of Cyg X-3. Also, the pairs of the pa- 
rameters (to, p) and (i, p) are each strongly anti-correlated. Still, 
our results indicate that, in spite of the strong X-ray orbital modu- 
lation, the inclination required to account for it does not need to be 
large, even for our spherically-symmetric wind model. The optical 
depths at the phases of 0, 0.25 and 0.5 are r = (1.5, 0.7, 0.5) and 
(1.0, 0.4, 0.3), respectively. Based on the fits, we can also calcu- 
late the corresponding mass loss rates using equation (|7), the def- 
inition of t , and the Cyg X-3 parameters adopted in Section |2~2l 
For the BAT and ASM best fits, r corresponds to -M = (0.8- 
1.4) x 10~ 5 [o"t/(ct c + <r ion )]MQ yr~' (where <x T is the Thomson 
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Figure 7. The fractional modulation depth, D(E), for all the available ASM 
and BAT data (blue and black symbols) and the ISGRI data (magenta and 
cyan symbols) estimated using fits with the model of equations (8), 111) 
(blue and magenta) and with the 3-harmonic model (black and cyan). The 
dashed curve shows the dependence of D(E) for equation (7} in the case of 
the Klein-Nishina cross section, approximately normalized to D(20keV). 



cross section), which is of the order of the values usually estimated 
for Cyg X-3 (see SZ08 and references therein). 

Fig-Oa-b) compares the folded light curves for the count rates 
of the ASM and BAT instruments for the hard, intermediate and 
soft states, using the criteria described in Section |2~T1 see also Fig. 
[2] We see that in the ASM range, the soft state modulation is sig- 
nificantly stronger than in the hard state, with the fractional modu- 
lation depth of D ^ 0.64, 0.57, respectively. A similar result for the 
ASM data (but using differ ent state criteria) is shown in fig. 6 of 
iHialmarsdotter et al.1 (12008). In the BAT range, the difference be- 
tween the states is smaller, D = 0.49 and 0.53 in the hard and soft 
state, respectively. The difference between the states may be due 
to the soft state taking place during periods of a higher wind mass 
loss rate. For the ASM and BAT, the ratio of the fitted values of To 
of our wind model between the soft and hard states is ^ 1.20 and 
= 1.08, respectively. The latter value may be close to the ratio of the 
mass loss rates (as the dominant opacity in the BAT range is likely 
to be Compton), whereas the former is likely to be also affected by 
changing the bound-free opacity due to the change of the ionizing 
continuum. 

We then compare the orbital modulation in different energy 
channels of the ASM and BAT, for all of the available data. We 
number the channels consecutively with the increasing energy, 
from 1 to 11. Fig. [6] shows some examples of comparison of the 
modulation profiles in different energy channels, and its panels are 
marked with the channel numbers being compared. We see that al- 
though the modulation depth generally decreases with the increas- 
ing energy, the strongest modulation occurs for the 3-5 keV range, 
as shown in Fig. [6ja). We have checked that the same behaviour 
takes place within the hard and soft states. This effect appears to be 
due to an additional, less or differently modulated, spectral com- 
ponent appearing in soft X-rays, which is probably re-emission in 
soft X-ray lines of the absorbed continuum by the stellar wind, 
see SZ08. The general decrease appears then to be due to the 
bound-free absorption cross section decreasing (except for ioniza- 
tion edges) and the Compton scattering cross section decreasing 
as well with the increasing energy. The profiles also become more 
symmetric with the increasing energy, indicating that the electron 
density is more symmetric (with respect to the plane perpendicular 
to the binary plane along the line joining the stars) than the ioniza- 




Figure 8. Comparison of the folded light curves in different energy chan- 
nels selected by spectral state. The blue, red and green light curves are for 
channels (a) 24-35 keV, 50-75 keV, 75-100 keV in the hard state, (b) 24-35 
keV, 50-75 keV, 100-150 keV in the soft state, and (c) 24-35 keV, 35-50 
keV, 50-75 keV during periods of high-energy y-ray emission, respectively. 



tion structure. All the profiles have the minima within <pj2n ^ 0.9- 
1, with no indication of a different modulation pattern up to 150 
keV. 

Fig. [7] shows the fractional modulation depth, D of equation 
d 1 2b . as a function of the photon energy, estimated using fits with 
the model of equations {SJ, i ll It with a free /3 and of the 3-harmonic 
model. Since the obtained statistical errors are very small, a mea- 
sure of the uncertainty is provided by the differences between the 
values for the two models. We show D only up to 75 keV, as at 
higher energies strong noise prevents trustworthy estimates. We see 
that indeed the orbital modulation is the strongest in the 3-5 keV 
range, and decreases towards both lower and higher energies. Fig. 
|7]also shows D from the ISGRI data (Section[3]below), which is in 
a good agreement with the BAT results. 

The dashed curve in Fig. [7] shows the dependence predicted 
using equation dl2t for Compton scattering alone, using the Klein- 
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Table 1. The log of the INTEGRAL ISGRI observations. The total exposure 
is about 4.5 X 10 6 s. The first row includes all the observations before the 
start of the BAT monitoring. The low 14-50 keV state corresponds to the 
soft state as defined in Section|2] 



Figure 9. Comparison of the folded light curves in the hard X-rays (red 
symbols) in the soft state in the 60-150 keV range with those at > 0.1 GeV 
(from FLC09, blue symbols). The constant background level (FLC09) has 
been subtracted from the LAT fluxes. The black curve shows the 100 keV 
modulation predict ed (Z12) using the anisotropic Compton model fitted by 
iDubus etalffcOlOl) to the GeV data. All the profiles are normalized by the 
respective average values. 



Nishina cross section and Att = 0.75, for which the model de- 
pendence crosses the observed one around 20 keV. We see that the 
decline of D above 5 keV cannot be accounted for by the decline 
of the Klein-Nishina cross section with energy, and photoionization 
absorption has to contribute at < 20 keV (in agreement with, e.g., 
SZ08). 

Figs. [8ja-b) show the modulation-profile dependence on en- 
ergy during the hard and soft state, respectively. We see that up to 
» 100-150 keV, at which energies the statistics become poor, the 
profiles have the minima around 0.9-1 phase, the same as at low 
energies and as for the entire data. An analogous figure for the in- 
termediate state looks similar and it is not shown. 

Fig. [He) shows the modulation for the data simultaneous with 
the observations of high-energy y-ray emission (see Section ETt . 
We see the same modulation pattern as for other data, and, in par- 
ticular, neither a change of the positions of the maximum and min- 
imum nor a flattening of the profile with respect to the profiles at 
other states up to at least 75 keV. At higher energies, we also see 
no indication of the presence of another modulation pattern, but the 
data become noisy. 

The blue symbols in Fig. [9] show the orbital modulation ob- 
served at E > 0.1 GeV by the Fermi LAT (FLC09). The black curve 
shows the modulation of the jet emi ssion predicted at 100 keV us- 
ing the anisotropic Compton model jDubus et al.l2010l) . from fig. 6 
in Z12. The 100 keV model is relatively similar to that at 0.2 GeV, 
though somewhat shifted to higher phases (Z12). These jet modu- 
lation profiles are shifted in phase with respect to those observed 
in X-rays, and have the maximum around the superior conjunction. 
We compare then the predicted modulation of the jet spectral com- 
ponent with that observed at 60-150 keV, red symbols. This profile 
has been obtained by averaging the BAT 100-150 keV profile with 
the INTEGRAL one for 60-100 keV (see Section[3]below), both in 
the soft state. Indeed, we see it is similar to the modulation pro- 
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files observed in other states and very different from either the jet 
contribution pred icted at 100 keV or observed at > 0.1 GeV. 

We note that lDubus et alj d2010h have determined the location 
of the y-ray emitting region at a distance ~ 2a along the jet orig- 
inating at the compact object. Then, X-ray emission from the jet 
would be much less absorbed by the wind than that originating at 
the compact object. Furthermore, if the jet X-ray emission were 
stronger in the soft state then the modulation amplitude would then 
decrease. Instead, we see in Fig. [8] that the modulation depth ei- 
ther in the soft state or for the data simultaneous with the observed 
y-ray emission is basically the same as in the hard state, during 
which there is no detectable jet contribution. We note that the mod- 
ulation pattern of the jet 100 keV emission can be flatter than that 
shown in Fig. [9] due to the electron cooling slower than at > 0.1 
GeV, causing the 100-keV emission region along the jet to be more 
extended than that in high-energy y-rays. Still, a major jet contri- 
bution would flatten the observed total modulation profiles, which 
is not seen. Thus, based on Figs.[8]{9](also Section[3]below), we see 
no indication of a contribution from the low-energy tail of the y-ray 
emission. That emission has to have a low-energy cutoff above the 
hard X-ray range. Indeed, the jet models with the minimum Lorentz 
factors of the accelerated electrons of 1300-1500, which are shown 
by the blue curves in figs. 5(a-b) in Z12, do satisfy this constraint 
and contribute very little at 100 keV. On the other hand, their mod- 
els with the minimum Lorentz factors of 300 and 700 predict too 
much flux at 100 keV and appear to be ruled out. 



3 ORBITAL MODULATION IN INTEGRAL DATA 

We have studied IN TEGRAL data from t he ISGRI detector (a part of 
the IBIS telescope, lUbertini et alll2003l) . We have used all the data 
public as of 2012 May 1, available at the INTEGRAL Science Data 
Centre (ISDC). We select only data with the off-axis angle < 5°, 
which corresponds to the fully coded part of the detector (the inner 
10° x 10° of the field of view). These data span MJD 52629-55542. 
INTEGRAL observations consist of pointings with a typical length 
of 2-3 ks. Our data set comprises of 1824 such pointings, for a total 
exposure of ^ 4.5 x 10 6 s. We consider only the 20-100 keV energy 
range since the data above 100 keV have relatively poor statistics. 
The data have been reduced using the Offline Scientific Analysis v. 
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Figure 10. Folded light curves for INTEGRALflSGRI observations, taking 
into account (a) all observations and (b) those corresponding to the soft 
state. The blue, red and green light curves are for energy ranges of 20—40 
keV, 40-60 keV, 60-100 keV, respectively. 



9.0 package provided by the ISDC dCourvoisier et al.ll2003l) . with 
the pipeline parameters set to the default values. The created ISGRI 
light curves have a resolution of 100 s. Table[T|gives the log of the 
observations. 

We consider first all the observations, and then only those cor- 
responding to the soft-state intervals with the low BAT flux, see 
Section |2TI We perform the analysis using the standard f tools 
software, in which the folded light curves are unweighted averages 
of the fluxes within individual bins. Our results are presented in Fig. 
[Tol We see no indication of a change of the profile shape towards 
high energies in either of the two cases considered, confirming our 
results in Section 1231 above. We have also calculated the modula- 
tion profiles for only the data simultaneous with y-ray detections. 
We have obtained results entirely consistent with those for the soft 
state, but with somewhat worse statistics. 

In Fig. [3t a )> we compare the 20-40 keV modulation profile 
for all the ISGRI data with the corresponding BAT profile. We see 
a very good agreement between them, which is consistent with the 
mostly overlapping energy coverage of the analyzed data. A 20- 
40 keV modulation pr ofile for an early observation of Cyg X-3 by 
ISGRI is presented in lVilhu et al.l ( 120031) , whose results are similar 
to ours but with substantially higher statistical fluctuations. That 
paper also presents the modulation profile from the monitoring by 
the CGRO/BATSE during 1991-2000, which is also of almost the 
same shape as those of our ISGRI and BAT profiles. 



4 RXTE PCA/HEXTE PHASE-RESOLVED SPECTRA 

Here, we study the RXTE Proportional Counter Array (PCA) and 
High Energy X-ray Transient Experiment (HEXTE) data sets used 



Table 2. The log of the RXTE observations. The exposures are for the sums 
of all the used observations within each phase interval within a state. 

Spectral state Inferior Conjunction Superior Conjunction 



Exposure [s] 



Exposure [s] 





PCA 


HEXTE 


PCA 


HEXTE 


1 


9120 


9120 


23408 


19549 


2 


24048 


14611 


42976 


25952 


3 


20912 


23049 


23008 


9337 


4 


3744 


4645 


2960 


3805 


5 


1120 


5465 


12688 


11672 



by SZM08, who obtained average X-ray sp ectra in five spectral 
states. The log of the observations is given in H ialmarsdotter et all 
(2009). We split each average data set into two subsets correspond- 
ing to the superior and inferior conjunctions, which is done by di- 
viding the orbital phases into two parts, tpjln = 0.28-0.78 and 
0.78-0.28. These boundaries have been determined us ing the or- 
bital template of Ivan der Klis & Bonnet-Bidaud 1 19891) with the 
condition that the orbital period is divided into two parts of equal 
length, with the template flux greater or lower than that correspond- 
ing to the dividing phases. 

We have performed standard RXTE PCA and HEXTE data re- 
duction using ftools. The resulting exposures for the two subsets 
of the data of SZM08 are given in Table[2] We see that whereas the 
states 1-3 have exposures > 10 4 s for each subset for either PCA 
or HEXTE, the states 4-5 have much worse coverage, especially 
for the PCA, with the exposures as short as a 1000 s and a 3000 
s. Therefore, we present only the spectral shapes for the states 4- 
5, whereas we study in detail the differences between the subsets 
around the superior and inferior conjunctions for the states 1-3. 

In our approach, we study the difference between the spectra 
around the inferior and superior conjunctions assuming it is due to 
an additional optical depth, with scattering and absorption charac- 
terized by a con stant ionization coefficient, see equation (B3) of 
IZdziarskil J2012l) . We first fit the spectra around the inferior con- 
junction. Finding a physical model of unabsorbed spectra of Cyg 
X-3 ha s proven, in this and p revious studies, to be difficult. For ex- 
ample, [HjdmMsdotteLetaL, ( 2008, 2009) found that fitting data at 
> 3 keV do not allow an unambiguous determination of the intrinsic 
spectra, with the range of possible models strongly varying at low 
energies. On the other hand, Cyg X-3 data from BeppoSAX extend 
to < 1 keV, and are better suited for such a determination. Those 
data have been studied by SZ08, who used a relatively complex 
model of wind absorption and emission , modifying spectra from 
hybrid plasma and Compton refle ction jPoutanen & Coppilll998t 
ICoppilll999HGierlinski et ai]|l999l) . 

Here, we fit the data from PCA and HEXTE, which are for 
E > 3 keV only, and do not allow unambiguous determination of 
the low-energy intrinsic spectra. However, our objective is to study 
physical mechanisms of the orbital modulation, and determining 
the actual intrinsic spectral distribution is beyond the scope of this 
work. We thus fit the data with a model that has been shown to pro- 
vide re latively good fits to a large set of data from RXTE a nd INTE- 
GRAL lVi\hu et alj2003l:lHialmarsdotter et al.l2008ll2009l SZM08, 
ICorbeletal .1120121) . but which still should be considered as phe- 
nomenological. In particular, any inferences about the unabsorbed 
spectra based on that model are very uncertain. The model is based 
on emission of hybrid plasma and reflection (including the Fe K 
line), similar to the model used by SZ08. However, the complex 
absorption by the wind is treated only phenomenologically, as a 
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Figure 11. The panels (a-c) show the spectra for the states 1-3, respectively, 
in the classification of SZM08. The upper (blue) and lower (red) spectra cor- 
respond to the phase ranges around the inferior and superior conjunction, 
respectively. The crosses show the unfolded PCA+HEXTE spectra. The in- 
ferior spectra are modelled by emission of hybrid plasma and reflection/Fe 
K line absorbed by a dual absorber. The models of the superior spectra are 
the best-fit models of the inferior spectra absorbed/scattered in an additional 
medium, corresponding to the difference in the average column density be- 
tween the two phase ranges. 



product of full and partial neutral absorbers. The Fe abundance of 
the absorber and reflector is assumed to be free. However, given the 
neglect of photoionization and the complex spatial structure of the 
wind, the fitted value should also be considered as phenomenolog- 
ical only. In spectral fitting, we allow a free normalization of each 
HEXTE spectrum with respect to the corresponding PCA one, and 




Figure 12. The spectra of the states 4-5, shown in the panels (a-b), respec- 
tively. The upper (blue) and lower (red) spectra correspond to the phases 
around the inferior conjunction and superior conjunction, respectively. The 
crosses represent the unfolded PCA+HEXTE spectra. The models (not 
shown) are given by emission of hybrid plasma and reflection/Fe K line 
absorbed by a dual absorber. 



we add a 1 per cent systematic error to each of the count spectra. 
The PCA and HEXTE data are fitted in the energy ranges of about 
3-25 keV and 20-40 keV, respectively (as we have found the sta- 
tistical quality of the HEXTE spectra above 40 keV is relatively 
poor). The resulting fits to the spectra in the states 1-3 are sta- 
tistically satisfactory, x 2 /v < 1- The fitted para meters are consistent 
with those pr esented in lVilhu et alJd2003h and Hialmar sdotter et akl 
d2008ll2009h . 

Then, we take the best-fit model to an inferior-conjunction 
spectrum, and fit the corresponding superior-conjunction one with 
the inferior-conjunction spectral model but attenuated by absorp- 
tion and scattering in an additional medium due to lines of sight 
through the wind around the superior conjunction being longer 
than those around the inferior conjunction. We have found that the 
fractional modulation depth averaged over our phase intervals is 
D = 0.4 but it strongly decreases at E>15 keV. We have first 
tested whether this behaviour can be due to the Klein-Nishina de- 
cline of the Compton cross section. However, that decline is much 
too slow to account for the observed energy-dependent modulation 
depth, which confirms our results in Section |2~31 above. Clearly, 
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Table 3. The fitted parameters of the additional absorbing medium related 
to the superior conjunction for the states 1-3 (which have the exposures 
sufficient for averaging). See Section|4]for details. 



State 




N pmt 


/a iv 


togio? 


x 2 /y 




10 23 cm 2 


10 23 cm 2 




erg cm s 




1 


2.0 ±0.2 


6.8 ±0.5 


0.36 ± 0.01 


1.6 ±0.4 


40/53 


2 


2.1 ± 1.4 


10.8 ± 8.7 


0.31 ±0.08 


3.3 ± 0.4 


30/55 


3 


3.5 ± 0.3 


3.0 ±0.6 


0.44 ± 0.07 


1.8 ±0.4 


55/56 



photoionization, known to be present in Cyg X-3 (e.g., SZ08), also 
contributes. Its cross section generally declines fast with increas- 
ing energy (apart from ionization edges). Accordingly, we have 
found that good fits to the data can be obtained including absorp- 
tion by an ionized medium. Specifically, our model for the addi- 
tional superior-conjunction attenuation consists of partial covering 
by a partially ionized medium and full covering by a fully ionized 
medium. Such a model may approximate shadowing of the stellar 
wind by the X-ray source, and it has been invoked t o explain the 
IR modulation and spectral lines dvan Kerkwiiklll993h . The attenu- 
ation is then by the factor, 

A(E) = exp(-/V fu ncrc) {l - / cov + / cov exp [-/V part (<x c + cr ion )]} , (13) 

where Nm and N pM are the electron column densities of the fully 
and partially covering medium, a nd /* cov is the coveri ng fraction. 
We use the ionized absorber of iReeves et alj d2008l) . based on 
xstar jBautista & Kallmar]|200ll) . It is characterized by the col- 
umn density and ionization parameter, £ = Lx/(w 2 ), where n is the 
medium dens ity and r is the distan ce from the ionizing source to the 
medium; see lReeves et alj d2008l) for further details. Our model is 
clearly very simplified. We have neglected spatial variations of the 
ionization coefficient. Also, we have assumed that scattered pho- 
tons are entirely removed from the line of sight, whereas some pho- 
tons are also scattered into the line of sight from other directions. 
The fit results are summarized in Table [3] and the unfolded spectra 
and the models are shown in Fig.fTTI 

We see that the total electron column density is /V ~ 10 24 cm 2 , 
which corresponds to the difference in the average Thomson optical 
depth between the two phase intervals of At t ~ 1 . This is in agree- 
ment with At ~ 1 found from fitting the orbital modulation profiles 
(Section[23}. The fitted values of £ ~ 10 2 ~ 3 erg cm s correspond 
to moderate ionization. We can notice that the models do not fit the 
smeared edges clearly seen around ~ 10 keV, which are apparently 
due to absorption by H-like Fe in a very strongly (but not fully) 
ionized phase of the wind, which we do not model here. Using the 
definition of f , the characteristic distance of the ionized medium is 
r ~ Lx/(Ng), which, for the fitted values of f and N is r ~ 10 11 " 12 
cm. This is of the order of the orbital separation, a, as expected. 
Thus, the fitted model is consistent with bound-free absorption tak- 
ing place in the stellar wind. We note that the fitted parameters do 
not form monotonic sequences with the increasing state number. 
This appears to be due to the limited observation exposures. Then, 
the two sub-states of a given state are not true averages, but are 
affected by variability taking place during the used observations. 

In Fig. [T2] we show the inferior and superior spectra for the 
states 4-5. We can see that due to the short exposures, the intrin- 
sic variability is apparently comparable to the orbital modulation, 
especially for the spectra of the state 4, which have the superior- 
phase fluxes at > 20 keV higher than those of the inferior phases, 
contrary to the expectations for long-term averages. 



5 CONCLUSIONS 

We have studied orbital modulation of X-rays from Cyg X-3. Us- 
ing the data from the detectors ASM, PCA and HEXTE of RXTE, 
BAT of Swift, and ISGRI of INTEGRAL, we have found that the 
modulation depth at E > 5 keV decreases with the photon energy. 
The decrease is too fast to be accounted for by the Klein-Nishina 
decrease of the Compton cross section with energy. As we confirm 
by fitting the PCA/HEXTE data, the energy dependence requires 
the presence of a moderately ionized absorber, with the ionization 
coefficient of LxKnr 2 ) ~ 10 2 ~ 3 erg cm s . We have also found 
a decrease of the depth of the modulation at E < 3 keV, which is 
probably due to the re-emission of the absorbed continuum in soft 
X-ray lines by the wind. 

Given that we use the data accumulated over long time, we 
determine the energy-dependent folded and averaged light curves 
with high accuracy. The modulation profiles are not consistent with 
a spherically symmetric wind. The measured detailed shapes of the 
orbital modulation and their photon-energy dependence can serve 
for determination of the wind structure in future work. Also, the 
energy dependence of the modulation depth at low energies can be 
further studied using spectra from detectors sensitive at soft X-rays. 

Assuming that the modulation is from wind absorption, the 
phase of the flux minimum corresponds to the maximum optical 
depth through the wind, which is then close to the superior conjunc- 
tion. On the other hand, high-energy y-rays detected from Cyg X-3 
by Fermi are found to have an orbital modulation shifted by about 
a half of the period with respect to that of the X-rays (FLC09). 
This pattern has been interpreted as due to anisotropy of scattering 
of stel lar photons by relativistic electrons in the jet by Dubus et al.1 
(2010). As calculated by Z12, the jet emission at hard X-rays has 
the modulation pattern almost the same as that in the GeV range. 
Thus, the jet X-ray contribution can potentially be measured by de- 
termining the orbital-modulation profiles at hard X-rays. 

We have performed such study, and have found no presence of 
an X-ray modulation pattern different from the standard one, i.e., 
with the minimum around the superior conjunction, up to ~100 
keV. This implies that an X-ray jet contribution is weak up to at 
least that energy and the jet component has to have a low-energy 
cutoff. This is consistent with the X-ray spectral measurements 
taken during periods of y-ray emission, which do no t show a hint 
of a j et component up to comparable energies (Z12: ICorbelet all 
|2012|) . This constraint is satisfied by two among the models pre- 
sented in Z12, those with the minimum Lorentz factors of the ac- 
celerated electrons of 1300 and 1500. 
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